Surface antiferromagnetism and incipient metal-insulator transition in strained manganite films 
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Using first-principles calculations, we show that the (001) surface of the ferromagnet Lai-^Sr^MnOa under 
epitaxial compression favors antiferromagnetic ordering in the first few surface layers, coexisting with ferromag- 
netic bulk order. The change in magnetic order is accompanied by a marked surface-related spectral pseudogap, 
signaling a metal-insulator near-transition at the surface. All surfaces exhibit huge outward relaxations an rum- 
pling of the MnC>2 and LaO surface planes, causing a distinct shift in work function. The AF phase is an extreme 
surface-assisted case of the combination of in-plane AF super-exchange and vertical FM double-exchange cou- 
plings that rules magnetism in manganites under in-plane compression. 

PACS numbers: 73.20.-r, 75.70.Ak, 75.47.Lx, 75.70.Rf 



Rare-earth manganites are promising compounds for inno- 
vative magnetic sensors (HE), spintronic devices |3-8|, and 
fuel-cell components l9l [T0ll . and have been the focus of many 
theoretical studies ITTTff4"ll . The colossal-magnetoresistant 
above-room-temperature ferromagnet (FM) Lai-^Sr^MnOa 
(LSMO henceforth) around £=1/3 is one of the most stud- 
ied in this class. Its magnetic phase and structural properties 
are quite sensitive to hydrostatic as well as uniaxial strains 
|[T3l [141 . It is natural to expect that such effects will be am- 
plified at surfaces and that tight magnetic phase competition 
may ensue, as suggested by previous theoretical work lfT31 - 
ITTl and by experiments on LSMO measuring a lowered mag- 
netic moment in substrate-constrained films |[T8l[T9l . Here we 
demonstrate, using ab initio simulations, that indeed antifer- 
romagnetism (AF) takes over at the (001) surface of LSMO 
under in-plane compression, and it coexists with FM order 
inside the bulk. The surface exhibits huge dipole-inducing 
rumplings which affect the work function, and its electronic 
properties suggest an incipient surface metal-insulator transi- 
tion and reduced in-plane conductivity. 

First-principles calculations are performed using the PAW 
method |20| and the GGA+U generalized gradient approx- 
imation with Hubbard correction |21| in the Dudarev et al. 
E2l formulation as implemented in the VASP code |23|. The 
only independent parameter, U-J, is set to 2 eV as in our previ- 
ous work on bulk LSMO lfT3l fl4l . GGA+U opens a gap at the 
Fermi level (Ep) in the minority channel [13, 14 24 1, in agree- 
ment with the known semimetallic character of LSMO. (In 
GGA, Ep barely slices through the lower conduction bands, 
so GGA and GGA+U results for structure and magnetism are 
in fact largely equivalent.) The energy cutoff is 400 eV, and 
Monkhorst-Pack k-points grids 4x4x1 and 9x9x2 are used 
for relaxation and density of states (DOS), respectively. 

We simulate the Mn-terminated (001) LSMO surface with 
a 13-layer symmetric slab with \^2x \[2 in-plane section, and 
containing seven Mn-02 and six (La/Sr)0 layers. The in- 
plane lattice constant is 3.80 A, i.e. that of LaAlOs, corre- 
sponding to an in-plane epitaxial compression of -2.6 %. We 
dope the material at £=1/3 by Sr substitution of one La per 
LaO layer, symmetrically with respect to the central layer; 



only the top LaO layers are Sr-free. The periodic slab im- 
age are separated by ~10 A of vacuum. The starting con- 
figuration is the bulk rhombohedral R3c with c/a= 1.05 (c is 
the (001) axis), but all ionic degrees of freedom are subse- 
quently relaxed. After an extensive search for competing sur- 
face magnetic structures, we restricted the comparison to the 
FM, q=(0,0,0), and to the so called "If" and "2f" structures, 
FM phases with Mn in the surface and, respectively, subsur- 
face Mn02 layer being spin-flipped according to an AFc pat- 
tern with q=( 1/2, 1/2,0) as in Fig[T| rightmost panel. Other 
patterns which are not energy-competitive are the 3f, i.e. the 
same as If or 2f with the first three Mn layers flipped, and 
other locally mimicking bulk AFa (q=(0,0,l/2)) and AFq 
(q=( 1/2, 1/2, 1/2)). 
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FIG. 1. (Color online) Simulation supercell (small red circles: O; 
large gray circles: La/Sr; intermediate blue circles: Mn). Left: top 
view (in-plane surface cell outlined). Center: side view (including 
vertical rumpling of surface layers). Right: If and 2f magnetic struc- 
tures (surface is upper side). 



From the raw total energies, we find that the If struc- 
ture, i.e. a planar G-type AF layer on top of the FM slab 
(Fig[T|, is favored over the FM. This shows unambiguously 
that under compressive epitaxial strain the first surface layer 
of LSMO (001 ) will be antiferromagnetic. The 2f phase is dis- 
favored, indicating that antiferromagnetism will hardly pene- 
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trate deeper than one Mn layer into the bulk. This is the case 
even in the best case scenario where we factor out the ener- 
getics of the bulk, using in lieu of the raw energy a surface 
specific |25| energy gain A7 = (AE S — nAEf,) /2A, where 
AE S is the energy difference between the FM slab and the 
spin-flipped slab, AEb is twice (there being two Mn per layer) 
the bulk energy difference per formula unit between AFc and 
FM, and n the number of layers which are spin-flipped in the 
slab (A=28.9 A 2 is the surface area). For the If structure, the 
gain is a sizable 1.8 meV/A 2 , while for the 2f structure the 
gain is zero to numerical accuracy (all other structures give 
negative gains). This confirms that only the topmost layer 
of the surface region of LSMO (001) will exhibit antiferro- 
magnetism at this strain level. Indeed, AF is expected to be 
surface-induced and surface-confined: previous analysis ||26l 
of bulk LSMO under strain has shown that the vertical mag- 
netic couplings are large and positive (i.e. FM) for all reason- 
able epitaxial conditions and external strains. A corollary is 
that larger strains will most probably not help AF to dive any 
deeper into the bulk. 
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FIG. 2. (Color online) Structural parameters vs position (top sur- 
face: left; bottom surface: right). Top: Mn-apical O interlayer dis- 
tance; center: octahedron rotation arouns s; bottom, Mn magnetic 
moments. 

The strained surface exhibits quite outstanding features for 
a metallic surface: i) a very large outward relaxation of the 
top layer (almost 20% compared to bulk interlayer spacing), 
associated with ii) an increased Mn magnetic moment at the 
surface, and iii) a large rumpling in the top MnC>2 and LaO 
layers. A selection of calculated quantities is in Fig|2] Mn-O 
octahedra behave nearly as rigid units, except in the surface 
and subsurface layers, where structural and magnetic behav- 
ior change in lockstep. The tilt angles of the octahedra off the 
surface normal s=(001) are essentially zero all over the slab, 
i.e. octahedra are perfectly aligned along s and laterally rigid. 
The large rotations of the octahedra around s are roughly con- 
stant in the bulk region, with a Mn-O-Mn angle a=160° as 
in bulk LSMO at the same strain; the rotation then increases 
in the subsurface layer (a=155°), and finally decreases in the 
surface layer (a=170°). Concurrently, the vertical (i.e. along 
s) distance between Mn in the top layer and O in the LaO layer 



just below increases to 2.2 A (+10% relative to the bulk 2.0 
A), while that between the subsurface Mn and O in the LaO 
layer below drops to 1.9 A (-5%). This correlates with the 
Mn magnetic moment, which drops from 3.4 /ig in the bulk 
to 3.3 /is in the subsurface layer (-6%), and then increases to 
3.7 fi B in the top layer (+9%). 

TABLE I. Rumpling (A) of the top surface layers. 
FM AF 5 
Mn0 2 0.235 0.170 -0.065 
LaO 0.265 0.255 -0.010 

In addition, we find that the Mn02 surface plane is rum- 
pled, as sketched in Fig[T]and summarized in Table|IJ with Mn 
atoms being strongly displaced upwards relative to O atoms. 
The same pattern, even in quantitative terms, appears in the 
first LaO layer from the surface. The rumpling is very large: 
the Mn-O (La-O) height difference is over two (over twenty) 
times the analogous Ti-O (Pb-O) displacement in the typical 
ferroelectric titanate PbTiOa 127*1 . This displacement gener- 
ates a large surface dipole pointing out of the solid into vac- 
uum. In particular, the increased near-surface rumpling at the 
FM surface relative to the AF surface (Table[lj effectively, the 
only significant structural difference between the two) is mir- 
rored in a reduced workfunction of the FM (W=4.47 eV for 
FM and W=4.67 eV for AF). Of course, this matches the usual 
trend of metal work function lowering with increased surface 
corrugation. Based on this result, a work function increase 
in a strained LSMO film would be a diagnostic evidence of a 
FM-to-AF magnetic phase change. 
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FIG. 3. (Color online) Local potentials and their periodicity-filtering 
running averages for FM (dashed) and If surfaces (solid) between 
surface (left) and slab center (right). Energy zero is the Fermi energy. 
Red (gray) line at the center is the difference <5V aV c of the If and FM 
averages (amplified by a factor 10). 

In further support of this argument Fig[3]displays the local 
potentials and their periodicity-filtering running averages for 
the FM (dashed) and If surfaces (solid). The local potential 
has troughs at the atomic positions (the deepest ones are Mn, 
the others La). The red (gray) line at the center is the differ- 
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ence <5V avc of the If and FM averages amplified by a factor 
10: clearly <5V avo is zero in the bulk, it starts being perturbed 
by the surface at the second Mn layer, it picks up decidedly 
from zero just below the top Mn layer as expected from the 
dipole difference, it grows steadily as the surface is crossed, 
and finally saturates to the work function difference (0.2 eV) 
in the vacuum. 

The electronic structure signature of this "ferroelastic" 
ifTTl distortion is a surface weakening of the directional d z 2- 
mediated bonding, as borne out by Figj4j discussed below. A 
similar rumpling of Mn02 planes (La-0 layer rumplings were 
not mentioned) at FM, unstrained LSMO (001) was reported 
earlier 1171 to decay gradually into the bulk. In the present 
case of in-plane compression this feature is much stronger and 
more surface-localized: our rumpling is twice as large (0.27 
A vs 0.14 A), and it occurs only in the outermost MnC>2 and 
LaO layers only, while it is an order of magnitude lower, i.e. 
negligible, in the other layers. 
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FIG. 4. (Color online) Orbitally resolved DOS for the FM (left) and 
If (right) magnetic order, at the surface (top) and in the innermost 
layer (bottom). Majority (minority) spin is positive (negative), the 
FM surface is fully metallic whereas the If has a much smaller DOS 
and in one channel only. The bulk DOS is identical in the two cases, 
demonstrating that the AF surface does not affect bulk magnetism 
beyond about 1 nm depth. 

Figj4] reports the orbitally-resolved DOS projected on Mn 
for the FM (left) and If (right) magnetic order, at the surface 
(top) and in the innermost layer (bottom). The bulk DOS in 
the lower panels is identical in the two cases, demonstrating 
further that the AF surface does not affect bulk magnetism be- 
yond at most 1 nm depth; also, half-metallicity (a gap in the 
minority channel) is preserved in all cases both in the bulk and 
at the surface. While the bulk is metallic with roughly equal 
participation of both e g orbitals, the FM surface is still fully 
metallic, but with mostly d x 2_ y 2 and only residual d z 2 char- 
acter at Ep- The standard DOS narrowing and upward shift 
at the surface (which would provide strong d z 2 metallicity) 
is prevented by the huge outward relaxation, which indeed 
is associated with the d z 2 peak dipping below E^. The AF 
If surface has a real gap in both the t2 g and d x 2_ y 2 chan- 



nels. On the other hand, the d z 2 weight moves below E^ 
as in the FM, but this demotion is incomplete and a small 
but appreciable DOS survives at Ep- We thus end up with 
a deep surface-related pseudogap, signaling an incipient but 
incomplete metal-insulator transition at the surface. This is 
still technologically relevant, since conduction through a con- 
tact on this surface in, e.g., a tunneling-magnetoresistance or 
field-effect device would clearly be affected by this 1-eV-wide 
pseudogap. 



0.5 


-0.5 
0.5 


-0.5 



1 1 1 1 1 1 


1 1 1 1 
/-. 
." \ 
I 

I 

1 

\ 


Mn 


* -y 

d , : 

L 




i 1 i 1 i 1 




o v ~" v ' 

a 

i.i.i, 


— p, 

-- Py 

1 


-6 -4 -2 2 4 
E - E p (eV) 



FIG. 5. (Color online) Orbitally resolved DOS for the If AF surface, 
projected on Mn in the top layer (top) and on apical oxygen O a in 
the subsurface layer (bottom). Majority (minority) spin is positive 
(negative). The residual metallic character is due to (non-bonding) 
overlap of Mn d z 2 with O a p z - The correspondence of other bonding 
peaks of Mn and O is also quite evident. 

The marginal (in the sense of being related to a pseudogap 
as in bulk Be, for example) character of surface metallicity is 
supported by the orbitally resolved DOS for one Mn in the top 
layer (Fig(5j top) and the corresponding apical oxygen O a in 
the subsurface LaO layer (Fig|5] bottom). The residual metal- 
lic character is clearly due to the (non-bonding) overlap of Mn 
d z 2 with O a p z . The correspondence of other peaks of Mn 
and O is also quite evident. Despite the absence of a proper 
gap, it is very likely that the conductivity will be strongly di- 
minished in this case due to the drop in E^ DOS. In closing, 
we note that the surface electronic structure just discussed is 
fully consistent with the strain effects on magnetic coupling 
mechanisms discussed in Ref 14 (see in particular Fig. 14, 
left panel): basically, in-plane compression favors rotations 
around the vertical axis s, and unbalances the magnetic cou- 
plings so that double-exchange (FM) tends to prevail along s 
while super-exchange (AF) is stronger in-plane (also, a charge 
transfer from d x 2_ y 2 to d z 2 attends to this shift). However, 
very large strains are needed in the bulk for in-plane AF to 
prevail. Here the surface provides locally the additional struc- 
tural freedom to tip the balance in AF's favor. 

In summary, we have shown that under compressive epi- 
taxial strain the first Mn surface layer (possibly, at most, the 
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first two) of ferromagnetic LSMO (001) is antiferromagnetic. 
This change in magnetic order is accompanied by a marked 
surface-related spectral pseudogap, signaling an incomplete 
metal-insulator transition at the surface. We found that the 
Mn02 and LaO surface planes are severely rumpled due to 
a weakened Mn d z 2-0 p z bonding. Despite metallicity, the 
rumpling-induced dipole causes a distinct work function shift. 
The surface AF phase is an extreme surface-assisted case of a 
previously recognized [14] combination of in-plane AF super- 
exchange and vertical FM double-exchange couplings under 
in-plane compression. 
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